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In this talk, we describe part of our recent work pQ (see also 12131 ) that gives new results in the context of 
asymptotic expansions of Feynman diagrams using the Mellin-Barnes representation. 



1. Introduction 

In the Standard Model quantum field theories, 
where calculations of Feynman diagrams are of- 
ten reduced to the evaluation of complicated in- 
tegrals involving several scale parameters, exact 
results are difficult to obtain. It is therefore legiti- 
mate to search to compute asymptotic expansions 
of these integrals in the various relevant limits 
of the parameters (or ratios of them) . A power- 
ful mathematical strategy to compute asymptotic 
expansions of integrals depending on parameters 
(in the limits of small and/or large values of the 
parameters) is based on the Mellin-Barnes repre- 
sentation |4I5) . The use of Mellin-Barnes repre- 
sentation in quantum field theory goes back in 
ref.[H] and many applications have been made 
since then. In ref. pQ we have proposed a new way 
for performing calculations within this strategy, 
using a theorem and new mathematical objects 
recently introduced in algorithmics 7 . In this 
talk, we shall briefly describe these mathemati- 
cal tools and show the power of our approach on 
a few diagrams of interest for the investigation 
of the anomalous magnetic moment of the muon. 
As a pedagogical example, a very elementary cal- 
culation will be exposed in a detailed way. 
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2. General lines of the strategy and some 
mathematical tools 

The underlying mathematical procedure of our 
method can be nicely summarized in only one sen- 
tence 7 : "For estimating asymptotically a func- 
tion f(x) [in the limit of small or large x], deter- 
mine its Mellin transform Ai[f](s) and translate 
back its singularities into asymptotic terms in the 
expansion of f{x)". 

Therefore, in the case of a Feynman diagram, 
we see that the problem is equivalent to get the 
Mellin-Barnes representation of its corresponding 
Feynman integral in a first step (since the Mellin- 
Barnes representation is the other denomination 
of the inverse Mellin transform) and once done, 
in a second step, to care for the singular structure 
of its explicitly obtained Mellin transform. 

The first step can be obtained in an easy man- 
ner: first use the Feynman paramctcrisation, then 
compute the momenta integrals to get purely 
parametric integrals and finally introduce an ap- 
propriate Mellin-Barnes representation of part of 
the integrand (or several if one has several con- 
venient parameters). Then, commute integrals 
over Feynman parameters and the Mellin-Barnes 
representation(s) and compute the parametric in- 
tegrals. In this way, the result is a Mellin-Barnes 
representation of the diagram where the Mellin 
transform of the diagram, directly readable in 
the integrand, is in general a product of Euler's 
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Gamma fonctions. In some of the examples that 
we have considered, the cosecant function also ap- 
peared. In these cases, it can be more convenient 
to keep this function explicitly rather than to ex- 
press it in terms of Gamma functions, since it can 
supply useful simplifications of the calculations. 

To deal with the second step of our problem, 
we need to define two important mathematical 
objects (the fundamental strip and the singular 
expansion of the Mellin transform of our diagram) 
and to give a theorem (the converse mapping the- 
orem). These tools have not been used previously 
in quantum field theory calculations but their ef- 
ficiency in this context was clearly shown in our 
work |1I2I3| . They are going to provide the dictio- 
nary between the singular structure of the Mellin 
transform of the diagram and its asymptotic ex- 
pansions. 

The fundamental strip < a, (3 > of the Mellin 
transform M [/] (s) of a locally integrable function 
f{x) (where x is a positive real number) is the 
largest open strip of the complex s-plane in which 
the integral 



M[f](s) 



dxx s 1 f(x) 



is convergent 

This fundamental strip can be easily found 
from the following conditions: 

If f{x) = 0{x~ a ) and fix) = 0{x~ p ) 

x— >0+ x— >+oo 

with a > (3 then M[f](s) converges in < a, (3 > 
(an alternative notation for the fundamental strip 
is »(a) €}a,/3[). 

The singular expansion of the Mellin transform 
in a given region of the complex s-plane, outside 
of the fundamental strip and where it admits a 
meromorphic continuation, is a formal sum col- 
lecting all singular elements of the Mellin trans- 
form in this region, a singular element (at a pole 
p) being simply the truncated Laurent expansion 
of the Mellin transform, at order O ((s — p)°) 
or less (depending on the Laurent expansion), 
around this pole 0. We note this formal series 



M[/]( S )x^a P) , 



(s-p) n 



Finally the converse mapping theorem is the 
following [7]. 

Let f(x) be continuous in ]0, +oo[ with Mellin 
transform A4[/](s) having a non empty funda- 
mental strip < a,/3 >. 

(i) If, for some 7 < a, -M[f](s) admits a mero- 
morphic continuation to the strip < 7, a >, if 
there exists 77 e]a, j3[ such that, in 7 < 5R(s) < r), 
M[f](s) = 0(\s\- r ) with r > 1, and if 

I s I — > 00 

Ai[f](s) admits the singular expansion 



M[f](s) i<J2 a p^ 
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in < 7, a >, then an asymptotic expansion of 
f(x) at is 
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+0(x-i). 



(ii) Similarly if, for some 7 > (3, M.[f](s) 
admits a meromorphic continuation to the strip 
< /3,7 >, if there exists 77 s]a, /3[ such that, in 
V < 3?(s) < 7, M[f]{s) = 0(\s\- r ) with 

s| — >oo 

r > 1, and if M[f}{s) admits the singular ex- 
pansion 



1 



(s~p') n 



in < /3, 7 >, then an asymptotic expansion of f(x) 
at +00 is 



p' ,n' 



(n' - 1)! 



x~ p In™ (x) 



+o( x - r ). 



In summary, each term of the left (resp. right) 2 
singular expansion of the Mellin transform of the 
function gives a corresponding term of the asymp- 
totic expansion of the function in the zero (resp. 
infinite) limit of x by simple reading. 

Some comments are in order. 
First, the fundamental strip of the Mellin trans- 
form is very important, because it separates two 
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distinct areas of the complex s-plane : a left 
and a right one. In the left region, the poles of 
the Mellin transform are linked, by the converse 
mapping theorem, to the asymptotic expansion 
of the original function f(x) in the zero limit of 
x. But poles in the right region are linked to 
the asymptotic expansion of f{x) in the infinite 
limit of x. It is therefore crucial to identify the 
fundamental strip of the Mellin transform in or- 
der to know which poles are associated to which 
asymptotic expansion. This is in general easy 
in the case of Feynman diagrams. As described 
above, when writing the Mellin-Barnes represen- 
tation of a diagram (after that the corresponding 
Feynman integral has already been transformed 
into purely parametric integrals), one uses the 
Mellin-Barnes representation of part of the para- 
metric integrand, so that one knows the funda- 
mental strip for this expression (using the rules 
explained above). Notice however that the evalu- 
ation of the parametric integrals, once their com- 
mutation with the Mellin-Barnes representation 
realized, can reduce the size of the fundamental 
strip. Therefore, the fundamental strip of the 
Mellin transform of a diagram is obtained only 
after evaluating the parametric integrals. 

Secondly, it is not necessary to compute the 
whole singular expansion(s) if one is only inter- 
ested in the first few terms, or in one particular 
term, of the asymptotic expansion(s). In these 
cases, one only has to consider the closest poles 
to the fundamental strip, or the pole correspond- 
ing to the particular term, in the relevant half 
complex plane. 

Finally, in many cases (in particular in all cases 
that we have considered so far), Mellin trans- 
forms admit a meromorphic continuation in the 
complete left and/or right half complex plane, 
therefore complete asymptotic expansion(s) for 
the original function can be found. Such expan- 
sions can be either convergent or divergent and, 
if convergent, they may in some cases represent 
the function exactly [7|. 

3. Applications: the anomalous magnetic 
moment of the muon 

The anomalous moment of the muon is a well- 
known precision observable for testing the Stan- 



dard Model (for a recent review, see jH])- 

One can easily use our method to compute the 
asymptotic expansions of the contributions to the 
anomaly of the diagrams of figure below (where 
t and t' are either electrons or taus with £ ^ £', 
and H denotes a hadronic polarisation 3 n#). We 
have therefore been able, among others, to di- 
rectly check in |1I2I3| the results of ^Oj (diagram 
2 of our figure) and eq. (8) of JT] corresponding 
to diagrams 4 and 5 (notice that eq. (8) of is 
not valid for a tau in the electromagnetic polarisa- 
tion). To be more precise, diagram 1 can be com- 
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puted in the —k — > +00 or —4- — > limits and di- 
agram 2 in both limits simultaneously. In the case 
of the other diagrams involving a hadronic polari- 
sation, perturbative and non-perturbative contri- 
butions are involved and it is therefore better to 
separate them with the use of a dispersion rela- 
tion written for the hadronic polarisation. One 
can then apply the asymptotic expansion proce- 
dure on the associated kernel function in the limit 
^2" - * +00 where t is the dispersive integration 
variable. 
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As a pedagogical example, let us consider one of 
the simplest of these diagrams in a detailed man- 
ner, namely diagram number 3. 

On can easily show that, once the dispersion 
relation has been worked out for II # , the contri- 
bution of this diagram to the anomaly is F2 (0) = 
a f^dt$(Il H (t))K 2 (^), where 

P x 2 m- 

is a well-known kernel function. We are interested 
in the asymptotic expansion of K% \ m?J wnen 

-T- — > +00. Although it is trivial to compute it 

"V 

from the exact result of the kernel function that 
is easy to obtain, let us show how to manage this 
evaluation within our strategy. 

3 Higher order hadronic polarisation contributions to the 
muon anomaly have first been studied in [5]. 
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First, introduce the Mellin-Barnes representa- 
tion 



l + X 



1 

2iTT 



7T 



c+zoo 

ds X~ 
c-too sin(Trs) 



(1) 



of fundamental strip < 0,1 >9 c and where X = 

x 2 'in' 2 , 



-\. One then gets 



A" 2 



2i7T 



c+ioo 



rf.S 



sin(7rs) 



x / dxx ia {l-xy- 8 . 
Jo 

This last equation involves the integral represen- 
tation of the Beta function B(2s + 1, 2 — s) which 
is equal to ^'ff+ff'' if 5ft(2s + 1) > and 
5R(2 — s) > 0. These conditions being always sat- 
isfied, since 5R(s) €]0, 1[ (the fundamental strip of 
eq. JU), we finally find the Mellin-Barnes repre- 
sentation of the kernel function 



finds 
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Kn 




with 



A n — 



1 T(2n + 1) 
(ra-2)! T(n + 3) 



B n = -A n [xjj(n - 1) - 2^(2n + 1) + ip(n + 3)] , 

where i/j is the Digamma function. 

Then, applying the converse mapping theorem, 
we finally have the asymptotic expansion 



K 2 



t 



1 

3~T 



( — 



^ Jc-ioo V m £ 



7T r(2s + i)r(2 - s) 

X sin(Trs) r(s + 3) ' 

Notice that in this case, the fundamental strip 
has not been perturbed after the integration over 
the Feynman parameter. The Mellin transform 
of the kernel function is then 



M 



t 



(s) = 



r(2s + i)r(2 



sin(Trs) T(s + 3) 



(2) 



of fundamental strip < 0, 1 >. 

Since we want the asymptotic expansion of K2 
in the infinite limit of -t-, we need to compute 

the singular expansion of @ to the right of the 
fundamental strip. Looking at J5J, we see that 
it can be easily obtained by calculating its trun- 
cated Laurent series around an arbitrary strictly 
positive natural integer n (the only singular func- 
tions of J2J in the right half complex plane being 
the cosecant and the T(2 — s) functions). One 



n=2 



A n 



In- 



{-l)B n 



Our strategy is, as we said, also useful for the 
other diagrams of the above figure. For diagrams 
2, 4 and 5, calculations are slightly more compli- 
cated but remain still easy to perform |1I2I3) . 
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